What Determines the Wave Function of Electron-Hole Pairs in Polariton 

Condensates? 
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The ground state of a microcavity polariton Bose-Einstein condensate is determined by consider- 
ing experimentally tunable parameters such as excitation density, detuning, and ultraviolet cutoff. 
During a change in the ground state of Bose-Einstein condensate from excitonic to photonic, which 
occurs as increasing the excitation density, the origin of the binding force of electron-hole pairs 
changes from Coulomb to photon-mediated interactions. The change in the origin gives rise to the 
strongly bound pairs with a small radius, like Frenkel excitons, in the photonic regime. The change 
in the ground state can be a crossover or a first-order transition, depending on the above-mentionsed 
parameters, and is outlined by a phase diagram. Our result provides valuable information that can 
be used to build theoretical models for each regime. 
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Microcavity polaritons — photocxcited electrons and 
holes strongly coupled with photons in a semiconduc- 
tor microcavity — have been observed to exhibit Bose- 
Einstein condensation (BEC) [l|, Due to the light- 
matter coupling, the polariton has an extremely small 
mass about 10^^ times the free-electron mass; the small 
mass results in a high critical temperature and low crit- 
ical density. BEC can be realized even at room tem- 
perature [3| , which is remarkable considering that it had 
been difficult to realize BEC in exciton systems for a long 
time [3l • Microcavity polaritons are dissipative particles 
due to the short lifetime of photons and inelastic scat- 
tering of excitons by phonons. Therefore, the polariton 
BEC is in a nonequilibrium stationary state with a bal- 
ance between pumping and losses 0,0, 13 ■ However, the 
polariton BEC has many similarities with BEC of neu- 
tral atoms in a thermal equilibrium Q . It shows several 
evidences for the superfluidity: the Goldstone mode 
the quantized vortices [l^l, and the collective fluid dy- 
namics fn\ . 

Such a stationary state appears to be well described by 
the ground state of a closed microcavity polariton system, 
when the polariton lifetime is longer than the thermaliza- 
tion time [6|,ll3| • In this paper, assuming such a situation, 
we determine the ground state with a fixed excitation 
density at absolute zero as a function of experimen tally 
variable parameters: excitation density, detuning [l3|, 
and ultraviolet cutoff determined by the lattice constant. 
In past studies, mean- field theories have been used to dis- 
cuss the two limits — low excitation density [3j and 
high excitation density — by considering two differ- 
ent models. These theories are complementary 



their relation is somewhat ambiguous. We investigate the 
intermediate density region as well, where the electron- 
hole (eh) wave function of the relative motion becomes 
important. We show that the ground state energy and 
wavefuction gradually change from those of excitons to 
photons as the excitation density increases. It is also 



shown that the change can be a crossover or a first-order 
transition, depending on the parameters considered. The 
latter transition is characterized by a jump in the pho- 
tonic fraction, and a sudden narrowing of the eh wave 
function that is accompanied by a change in the binding 
force from Coulombic force to photon-mediated force. 

The polariton system can be described as electrons and 
photons with a total excitation number iVex, which are 
interacting through electric dipole coupling. The Hamil- 
tonian is given hy H = + i?ci-ci + ^^ph + H^-ph [H, 
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Uq [= {2i:e? / e*Vq'^)\ are the electronic dispersion 
in an effective mass approximation and the Coulomb 
interaction, respectively. Further, a^, 5fc, and il^k 
are the annihilation operators of the conduction and 
valence electrons and photons with momentum fc, 
respectively. Fourier transform of the density operator 
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bk-qb\). The zero-point 



frequency of the cavity photons is ojc and detuning is 
defined as d = {hwc — Eg)/so, where Eq is the exciton 
Rydberg. The light-matter coupling constant is given 
17| . but hy g = dcv y^2TrhLLic/e*V. The momentum dependence 
of the dipole coupling is neglected here. Instead, a 
momentum cutoff kc is introduced so as to restrict the 
electronic states contributing to the polariton formation 
to |fc| < fee. It is smaller than or roughly equal to the 
inverse of lattice spacing (e.g., 60/ao for a GaAs-based 
microcavity, with uq being the exciton Bohr radius). 
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FIG. 1: Energy per excitation is plotted as a function of 
J?a: (a) d = 8,5,3,1,0,-0.7 with fccfflo = 20, g = 0.1 (solid) 
and (b) fccOo = 20, 100, 200, 300, 400 with d = 3 and § = 
0.1 (solid). The dashed curves in both figures are obtained for 
eh systems [3 • The dotted curve in (a) denotes the Hartree- 
Fock energy for eh plasma. 



Considering the coherent state of polarizations and 
photons, the mean-field ground state of a polariton 
condensate is given by 

1$) = e^^^o-^^"'>l[{e'^'uk + Vkalb,)\vSic), (5) 



where |vac) denotes the vacuum state with no conduc- 
tion electrons, no valence holes, and no excited pho- 
tons. There is a normalization condition, u\ -\- = 1, 
that allows us to include the phase-space filling effects of 
fermions. The variational parameters A, Xfc, Ufc, and Vk 
are determined by the minimization of the total energy 
E {= {H + /iiVcx)) for a fixed A^cx, which is given by the 
expression X]fc(''/'l;V'fc + (afcOfc + ^fc^fe)/2). In the coherent 
state, all the eh pairs are found to have the same phase: 
Xk = [31 • After angular integration, the mean-field en- 
ergy per excitation e(= E/Ncx) and the total excitation 
density Pox(= A^ox/^) are given by 
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where n = fcag, Kc 



(7) 



is the mean separation, A is the normalized photon field 

given by A = X^/Aiw^Jw , and g g^/wJAm]^epj is 

a dimensionless coupling constant. It is difficult to de- 
termine an infinite number of variational parameters, we 
apply for the excitonic constituent the inter pol ating wave 
function proposed by Comte and Nozieres |18j : 
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FIG. 2: Photonic fraction plotted as a function of Rb for the 
same parameter sets as in Figs.[T][a) and[ljb), respectively. 



Our task is to determine three parameters C, f^, and A. 

Figure [1] shows the mean-field energy per excitation 
plotted as a function of Rs- In Fig.[Tl^a), each curve (col- 
ored lines) corresponds to a different value of d; however, 
all the curves have the same cutoff parameter fccOo — 20. 
The curves approach that of an eh system (dashed) in 
the low-density region, except near resonance, where 
d ~ — 1 {hhJc Eg — £{)), and in the high-density re- 
gion, slightly below the photon level d. The energy shift 
from e/eo = —1 to e/eg = d indicates the polariton con- 
densation crossovers from excitonic to photonic ones with 
an increase in the density. The energy saturation at high 
density is explained by fermionic phase-space filling. The 
number of eh pairs increases until the conduction electron 
band is filled up to the photon level, and thereafter, pho- 
tonic excitations replace those of eh pairs to minimize the 
total energy. The photonic character is observed above a 
density (i?s < Rt) where the curves move away from the 
dashed curve. The crossover can be seen directly from 
the plots in Fig. [2l Figure [2ja) is obtained for the same 
parameter set as Fig. [Ija). The photonic fraction in- 
creases sharply from zero for Rs < R*, except for d ~ — 1 , 
where a large photonic fraction is present at low densi- 
ties. Slopes at Rs < R* are higher for larger values of d. 
The evolution of the state from excitonic to photonic is 
always smooth for the cutoff value of kcUo = 20. 

In the excitonic regime where Rg > Rs, the ground 
state is classified into three kinds even though there are 
no clear boundaries. If one-excitation energy of eh sys- 
tems (dashed in Figs. [TJa) and (b)) is close to the sigle- 
exciton level — eo, the eh pairs can be regarded as the 
BEC of excitons. The low-density regime {Rs > 2) 
is categorized as "exciton BEC." The highest density 
regime {Rs < 1) is categorized as "eh plasma" [31 since 
the curve of the eh system (dashed curve in Figs. [TJa) 
and (b)) overlaps with that of eh plasma (dotted curve 
in Fig. [Ha)), where eh pairs are unbound and all ex- 
citations are fermionic. For the intermediate-density 
(1 < i?s ^ 2), which is categorized as "eh BCS," eh 



pairs are regarded as weakly bound fcrmions like Cooper 
pairs of BCS superconductivity. Since R* varies, d deter- 
mines the regimes the polariton system passes through — 
exciton BEC, eh BCS, and eh plasma — before the ground 
state of the system becomes photon-like. 

The plots shown in Figs, [ija) and HUa) arc shown in 
Figs. [Ijb) and[2ljb) for a different paramter set i.e., for 
d = 3 and several values of kc. The change in the polari- 
ton state at /ccOo = 100,200,300 is not a crossover but 
rather a first-order transition since there is a discontinu- 
ous jump in the slope of the curve in Fig.[Tl^b). Similarly, 
as shown in Fig.[2llb), the photonic fraction jumps to al- 
most 100% at the transition density. The system ground 
state is described by another type of solution different 
from the one that approaches the ground state of the eh 
system at low densities. We shall call the former photon 
solution since the photonic fraction is almost 100%, and 
the latter polariton solution. For kcOo = 400, the photon 
solution corresponds to the ground state for all densities. 

As discussed above, the evolution of the ground state, 
with increasing excitations, depends on d and kc- Vari- 
ous types of the evolution are summarized in the phase 
diagram of Fig. [3l The labels X, x, eh, and pol denote 
the different regimes where the polariton solution corre- 
sponds to the ground state: exciton BEC, eh BCS, eh 
plasma, and polariton BEC, respectively [20|. The la- 
bel ph stands for the regime where the ground states are 
photon-like. Therefore, for the parameter regime marked 
as "X-x-pol/ph" at the center of Fig.[3l the condensation 
changes from/to exciton BEC, eh BCS, polariton BEC, 
and photonic BEC in increasing order of excitations. The 
hyphen (-) and slash (/) indicate that the change between 
the different regimes is a crossover and a first-order tran- 
sition, respectively. Boundaries between different phases 
are determined under certain conditions [20| . The bound- 
aries that depend on the conditions are shown as dashed 
lines. The diagram indicates the following: (i) The larger 
the value of kc, the stronger the photonic nature of the 
condensate, (ii) The lesser the value of kc and d, the 
stronger the coupling between the eh system and pho- 
tons, i.e., the more polaritonic is the condensate. 

To scrutinize the difference between the polariton solu- 
tion and the photon solution, we plot in Fig. [4] the wave 
function of an eh pair, P{r) = {l/V) J^ki'^k^k) exp(ifcr), 
as a function of the relative coordinate r between an elec- 
tron and a hole. Plots are obtained for various densities 
(J?,) for d = 3 and kcQo = 20 [Fig. IHa)] and for d = 3 
and kcQo — 100 [Fig. IHb)]. In the former figure, the 
wave function changes its shape and narrows gradually 
with increasing density from the wave funtion of Is exci- 
ton, indicating that the parameter used ( in Fig. [3]) is in 
the regime where all the expected changes in the ground 
state are crossovers. While the width oq is determined by 
the Coulomb attraction in the exciton BEC regime, it is 
modified by the electric dipole interaction for Rs < R^, 
where the photonic character appears. The change in 
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FIG. 3: Phase diagram showing the crossover properties of 
the polariton BEC when the excitation density increases. The 
coupling parameter is set as g — 0.1. Phase boundaries that 
are clear and unclear are shown by the solid and the dashed 
lines, respectively (see the text for an explanation). 



the binding force is seen more clearly in Fig. Hfb) ( in 
Fig. [3]). The wave function gradually narrows as Rg de- 
creases from 3 to 0.8, and it shows a drastic change to 
have a sharp peak at r = for Rs < 0.7, i.e., in the pho- 
tonic BEC regime. Clearly, the width is determined by 
a new length scale and not by oq, which indicates that 
the mechanism of eh pairing is completely different from 
that in the case of dilute excitons. We find that the pho- 
ton solution satisfies the conditions Ufe <C 1, Ufc ^ 1, and 
A 7^ 0. Therefore, the variational equation for the ch 
wave function P{r) reduces to 
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P{r) ~ g\~S{r) = ^,Pir), (9) 



where (5(r) — {^/V)J2\k\<k cxp(ifcr) is localized at ?- = 
with a width of ~ l/^c and becomes a delta function 
if /cc = oo- The reduced mass of an eh pair is denoted by 
rrir. Since the third term on the left side can be rewritten 
as — ^^^^_^^ (5(r)P(r), we see that the photon field induces 
an attractive delta potential mediated by photons [l6| . 
When considering two or three dimensions, the lowest 
bound-state energy of the real delta potential becomes 
minus infinity. However, a momentum cutoff was intro- 
duced here, hence the bound-state energy remains finite. 
The energy is estimated as e/eo ^ d — ^g^fccOo, which 
is obtained using some approximations and the condi- 
tions ffc <C 1, ^ — oo, and i?s^^|C/^l ^ 1 (the photon 
solution satisfies them). We can conclude that the first- 
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cavity polariton system is determined by a variational 
approach The ground state changes from excitonic 
to photonic regime. In the photonic regime, eh pairs 
are shown to have a small radius since they are bound 
by photon-mediated delta attraction. With increasing 
excitations, various types of ground-state evolutions are 
expected depending on the detuning and momentum cut- 
off. The phase diagram in Fig. [3] summarizes the evolu- 
tions; the change is a crossover when the cutoff momen- 
tum is not too large, and the first-order transition from 
excitonic to photonic condensation is expected when the 
cutoff is large. Although the calculation presented here 
is for three dimensions, the same scenario applies to two 
dimensions, which will be shown elsewhere. 

The authors are grateful to Kenichi Asano, Takuma 
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FIG. 4: Real space profiles of the eh wave function 
P{r)/R^'^^^ for (a) = 3,2,0.8,0.5,0.3 with fc.ao = 20 and 
(b) i?s = 3,2,0.8,0.5 with fccfio = 100. The inset shows a 
magnified view of the region r/ao < 0.5. We set g = 0.1 and 
d = 3 for both graphs. 



order transition occurs when the bound-state energy due 
to the short-range attraction falls below the energy of the 
polariton solution. The solution of Eq. ^ for large fcc is 
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where ip^^^ and E^^ ((^^ and E^) are the wave functions 
and the energy of ns-exciton bound state (scattering state 
with momentum k). As discussed above, the chemical 
potential has an upper limit, i.e., ^/e < d. For a large 
momentum, the main contribution to P{r) comes from 
the scattering states with P{k) oc fc^^, and P{k) is long- 
tailed in the photonic regime compared to the case of the 
dilute limit of Is excitons [P{k) oc k~'^]. This directly 
leads to the narrowing of P{r). 

Our result can give clues to build thcrctical models 
that can appropriately describe the BEG. Since the mo- 
mentum cutoff introduced here will be of the same order 
as the inverse of the lattice constant, the width of the eh 
wave function can be of the same order as the lattice con- 
stant in the photonic regime. In such a case, an eh pair 
should be recognized as a Frenkel exciton, and hence, the 
description is beyond the capability of our model which 
employs the effective mass approximation. This indicates 
the need for other theoretical models to treat localized 
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